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We study a system of ultra-cold fermionic polar molecules in a two-dimensional square lattice 
interacting via both the long-ranged dipole-dipole interaction and a short-ranged on-site attractive 
interaction. Singlet superfluid, charge density wave, and supersolid phases are found to exist in 
the system. We map out the zero temperature phase diagram and find that the supersolid phase 
is considerably stabilized by the dipole-dipole interaction and thus can exist over a large region of 
filling factors. We study the melting of the supersolid phase with increasing temperature, map out a 
finite temperature phase diagram of the system at fixed filling, and determine the parameter region 
where the supersolid phase can possibly be observed in experiments. 
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I. INTRODUCTION 

Ultracold atom physics has undergone spectacular de- 
velopment in the last two decades. During the last years, 
after the first theoretical proposal to simulate the Bose- 
Hubbard model with ultracold bosons in an optical lattice 
[l[ and its successful experimental realization Q, con- 
siderable theoretical and experimental progress has been 
made on degenerate quantum gases in optical lattices Q . 
One major common goal shared by these efforts is to sim- 
ulate typical systems or minimal models of condensed 
matter physics, e.g. the Hubbard model, using ultracold 
atoms confined in optical lattices (1, It is expected 
that these investigations can give key insights into unre- 
solved open questions in condensed matter physics, e.g. 
the mechanism of high-T c superconductivity. Usually, 
interactions in cold gases are isotropic and short ranged. 
Naturally one can ask the question which type of new 
physics would arise in these systems if the atoms have 
additional long-range anisotropic interactions. 

As a matter of fact, a wealth of interesting physics 
has already been revealed by theoretical and experimen- 
tal investigations of ultracold dipolar Bose gases both 
in harmonic traps [ol-Hij and optical lattices |15l - [l8| . In 
contrast to dipolar Bose gases, due to the Pauli exclu- 
sion principle, dipolar effects in fermionic systems can 
manifest themselves only if the dipolar interaction en- 
ergy is at least of the order of the Fermi energy. Although 
there have been already some theoretical investigations 
on these systems |19l - [23| . this condition could hardly be 
reached experimentally until the recent successful realiza- 
tion of a degenerate quantum gas of fermionic 40 K 87 Rb 
polar molecules [24]. With this achievement, the door to- 
wards exploring the many-body physics originating from 
dipole-dipole interactions in fermionic systems has been 
opened, which has already lead to a number of theoreti- 
cal investigations on these systems [25-43]. Motivated by 
experimental progress on ultracold polar molecules, here 
we investigate the low temperature quantum phases of 
a system of cold polar fermions loaded into a 2D square 
optical lattice. 

The rest part of this paper is organized as follows. In 



Sec. II, we introduce the system and the model studied 
here. In Sec. Ill, we describe the theoretical approach 
used in our investigation. In Sec. IV, the main section of 
this paper, we present a detailed discussion of the quan- 
tum phases at both zero and finite temperature. Finally, 
we conclude in Sec. V. 



II. SYSTEM AND MODEL 

Motivated by ongoing experiments, we focus our inves- 
tigations on possible setups based on 40 K 87 Rb molecules 
which have a permanent electric dipole moment of 0.57 
Debye (D), where ID = 3.336 x 1(T 30 C • m. However 
the effective dipole moment in the laboratory frame is 
zero in the absence of an external electric field. When an 
external electric field is applied the molecules align with 
the field and have an induced dipole moment d which in- 
creases with the strength of the external field. Currently, 
the experimentally accessible range of d is [0, 0.22 D] |44| . 
In this work we assume the external electric field is ori- 
ented perpendicular to the optical lattice plane. 

Loading the molecules into an optical lattice, the 
physics of this system for sufficiently low filling can be 
captured by an extended Hubbard model within the low- 
est band approximation, 

H = -t ]T clc^ + Uj^n^nn (I) 

<i,j>,a i 



n,-, 



Here, denotes nearest-neighbour sites, c| CT (ci CT ) is the 
creation (annihilation) operator of a fermionic molecule 
with spin er on site i in the Wannier representation, and 
rii a = c\ a Ci a is the particle number operator. The first 
term in Eq. JTJ) describes the kinetic energy with t be- 
ing the hopping amplitude; the second term represents 
the onsite interaction U between molecules with oppo- 
site spin; the third term originates from the long-range 
dipole-dipole interactions, where c c i = d 2 / (Airsoa 3 ) char- 
acterizes the dipole-dipole interaction strength with a 
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and £o being the lattice constant and the vacuum permit- 
tivity respectively, and Rj is the dimensionless position 
vector of lattice site i; finally, in the last term /i denotes 
the chemical potential. 

In our calculations we assume a 2D optical lattice cre- 
ated by laser beams with wavelength A — 1064nm as used 
in the experimental setup of Ref. [44[ . We consider the 
following parameters: the height of the lattice potential 
in the direction perpendicular to the lattice plane V\ &t 
and parallel to the plane V[ at are chosen as V_} at = lOEn, 
Vl at = 10U[ at , where E R = h 2 /(2m\ 2 ) is the recoil en- 
ergy, m is the mass of molecule and h is Planck's con- 
stant. Under these conditions, the dipolar interaction 
strength is in the range [0, O.OSEr], and the hopping am- 
plitude is approximately 0.2Er. Therefore the ratio be- 
tween the dipolar interaction strength and the hopping 
amplitude Cd/t is approximately in the range [0,0.15], 
which we assume in our calculations. 

A good estimate of the strength of the onsite inter- 
action U requires a detailed study of all types of short- 
range interactions between the polar molecules, which is 
beyond the scope of this work. Since the main aim of 
this work is to investigate the possibility of observing a 
supersolid phase of polar molecules, we assume that the 
onsite interaction is attractive, i.e., U < 0. Moreover, we 
notice that since large attractive on-site interactions will 
make the system unstable, the region of |J7| investigated 
in this work is restricted to [t, 8t]. 

On the route towards observing the supersolid phase 
in experiments, the temperature plays a dominant role. 
In current experiments with polar molecules in a har- 
monic trap, the lowest temperature which can be reached 
is of the order of the Fermi temperature Tp |44j|. Also, 
when the lattice potential is ramped up, generally, the 
molecules will be heated in this process, which results 
in even higher temperatures. In the following, the sys- 
tem is studied at both zero and finite temperatures, and 
the temperature region where a supersolid phase can be 
possibly observed is determined. 



III. METHOD 



we notice that in the high-dimensional limit inter-site in- 
teractions only contribute on the Hartree level [43 |. In 
other words, the Hartree term of the inter-site interac- 
tion will dominate as the spatial dimension of the system 
increases. This motivates us to keep only the Hartree 
contribution of the inter-site interaction in the Hamil- 
tonian as an approximation to the original Hamiltonian 
©, i.e. 



~2 , |R,-R 3 f 



(2) 
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Our DMFT Hamiltonian therefore takes the form 
H = -t ^ C L C 3° + U'^2n it n t i + '^2(V i - \i)n ia (4) 



where 



V = 
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and we have rescaled the chemical potential /i in Q ac- 
cording to the trivial constant term in (J3]). 

We investigate this system using real-space DMFT (R- 
DMFT)[H, Hi, which is an extension of DMFT to a 
position-dependent self-energy and fully captures the in- 
homogeneity of the system. Within DMFT/R-DMFT, 
the physics on each lattice site is determined from a local 
effective action which can be captured by an effective An- 
derson impurity model |47|. In this work, we use Exact 
Diagonalization (ED) [50L |51| of the Anderson impurity 
Hamiltonian to solve the local action. Details of the R- 
DMFT method can be found in previous works [IH |4(| . 

Within each R-DMFT iteration, the Hartree contribu- 
tions Vi are calculated and the Hamiltonian (j4]) is up- 
dated according to the new values V,. This modified R- 
DMFT iteration is repeated until a convergent solution 
is obtained. 



Before discussing our results in detail, we give a brief 
description of our theoretical approach. In polar molecule 
systems, we have to consider not only the on-site attrac- 
tion of the molecules but also the inter-site repulsion due 
to the dipolar interaction. For treating the on-site in- 
teractions, the Dynamical Mean-Field Theory (DMFT) 
[47l | is well suited since it is non-perturbative, captures 
the local quantum fluctuation exactly, and gives exact 
results in the infinite-dimensional limit. For long-ranged 
interactions in fermionic system, on the other hand, there 
are, to our knowledge, up to now only few efficient ways 
to treat them, although DMFT has already been gener- 
alized to the so called cellular DMFT in order to treat 
short-range inter-site interactions (48| . 

We include long-range interactions as follows. First 



IV. RESULTS 
A. Zero temperature 

1. Phase diagram 

At zero temperature, the main properties of the sys- 
tem are summarized in the phase diagram in Fig. [1] 
which shows the phase boundary between the supersolid 
(SS) phase and the homogeneous singlet superfluid phase 
(HSF) at different onsite interaction strengths U in terms 
of the filling factor p = J2ia( n icr) /N and the dipolar in- 
teraction strength q, where N is the number of lattice 
sites. This phase diagram is based on calculations for a 
12x12 square lattice. We will denote the phase boundary 
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between the supersolid (SS) phase and the homogeneous 
singlet superfluid phase the SS-HSF boundary in the fol- 
lowing. 

Along the vertical line at half filling (p — 1) in Fig. [TJ 
irrespective of the strengths of on-site and inter-site inter- 
actions, the system is always in an incompressible charge 
density wave (CDW) phase, which is characterized both 
by vanishing compressibility dp/ dp, = and a finite den- 
sity modulation in real space. More specifically, the den- 
sity distribution over the lattice has a checkerboard (CB) 
structure, which is characterized by pq =^2j e ' * Pj/N 
with pj = J2a( n j<?) an d Q = l 71 ": 71 ")- I n our calculations, 
we choose pq as the CDW order parameter and we will 
denote this incompressible CDW phase as "CB solid". 

Away from half filling but still in the region bounded 
by the SS-HSF boundary and the vertical line of half fill- 
ing, a supersolid phase is obtained by doping the CB solid 
with either particles or vacancies, which is characterized 
by the coexistence of the singlet pairing order parameter 
Aj = (ciic.if) and the CDW order pq. As we can see 
from Fig. [TJ at fixed on-site interaction strength U, the 
filling region in which the system remains supersolid in- 
creases with the strength of the dipolar interaction. This 
can be understood by the following simple argument. 
We notice that near the SS-HSF boundary the dominant 
contribution to the inter-site interaction energy between 
two nearest neighbor sites i, j is positive and approxi- 
mately proportional to CdPiPj ~ Cd{p 2 — (pi — /9j) 2 /4), 
therefore a larger a will make the system more likely 
to favor a modulated density distribution since it can 
lower the intersite interaction energy in this way. On 
the other hand, at fixed dipolar interaction strength 
Crf, we can see from Fig. [1] that a larger on-site attrac- 
tion strength |?7| corresponds to a larger supersolid fill- 
ing factor region. This can be understood by a simi- 
lar argument as above. We notice that near the SS- 
HSF boundary the total on-site interaction energy of 
two neighboring sites is approximately proportional to 
-PM + p]) ~ -\U\{p 2 + ( Pi - p J ) 2 /4), one can thus 
easily see that for large \U\ the system will favor a large 
density imbalance between neighboring sites which re- 
sults in a modulated density distribution. 

Concerning the proposed experimental setup with 
40 K 87 Rb molecules in an optical lattice, according to the 
discussions above, the supersolid phase will be observed 
more easily in an interaction region with relatively large 
dipolar interaction and onsite attraction, since this will 
stabilize the supersolid in a large filling factor region, 
e.g. in the range p € (1.0,1.4) for interaction strengths 
c d = O.li and U = -8i. 

On the right-hand side of the SS-HSF boundary, the 
system is in a homogeneous singlet superfluid phase 
which is characterized by a uniform distribution of both 
the density pi and the singlet paring order parameter A j . 

We note that in the classical limit of zero hopping, the 
long-range interacting model exhibits a devil's staircase 
of various solid phases (52| . Therefore, for much larger Cd, 
which far exceeds the region investigated in this work, one 
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Figure 1: (Color online). Phase diagram of the system with 
respect to the filling per lattice site p and the strength of 
dipolar interaction a, obtained from calculations on a 12 x 12 
square lattice. Different curves indicate the phase boundaries 
between supersolid (SS) and homogeneous singlet superfluid 
(HSF) at different onsite interaction strengths U. On the 
left-hand side of each curve the system is in the SS phase 
(except exactly at half filling (p = 1.0)) while the right-hand 
side region indicates the homogeneous singlet superconduct- 
ing phase. From left to right, different curves correspond to 
U/t = -1.0, -2.0, -4.0, and -8.0 respectively. 



would expect to find other types of incompressible CDW 
phases rather than CB solid and possibly new SS phases 
similar to those found for dipolar Bose gases in optical 
lattices [l7|, [HJ • A detailed investigation of this "devil's 
staircase" in the quantum case is beyond the scope of this 
current work. 



2. Density and sinylet pairing order parameter distribution 

Recently, there has been considerable progress in 
single-site addressability in optical lattices using elec- 
tron and optical microscopy which allows for a direct, 
in situ, experimental observation of particle positions 
and density-density correlations of the system (53T-|56|. 
In Fig. [5] we show three sets of snapshots of the density 
and pairing order parameter distributions of the system 
for the different phases discussed previously. We observe 
that the system develops the (ir, tt) density modulation 
in both the CB solid and SS phase. Moreover in the SS 
phase the amplitude of the pairing order parameter Aj 
also has a similar spatial checkerboard modulation. In 
the HSF phase both the density and the pairing order 
parameter are constant over the lattice. We remark here 
that although the distribution of density and pairing or- 
der parameter is obtained in the absence of the external 
trapping potential, the typical feature of these phases can 
nevertheless be observed at the center of a shallow har- 
monic trap where a local density approximation is valid. 





Figure 2: (Color online). Zero temperature density pi (left 
panel) and singlet pairing order parameter A, (right panel) 
distribution of the system (U = — 4.0t, Cd = O.li). In descend- 
ing order, the filling factor has the values p = 1.0, 1.28, and 
1.46 respectively, corresponding to the checkerboard solid, su- 
persolid, and homogeneous singlet superfluid. 



B. Finite temperature 

In this subsection we investigate finite temperature ef- 
fects and calculate the critical temperature of the super- 
solid which could be useful as a guide for future experi- 
ments. 

First, we study the melting process of the supersolid to 
a normal phase with increasing T at fixed filling. Fig. [3] 
shows the temperature dependence of the CDW order 
parameter pq and the average value of singlet pairing 
order parameter A = Aj/N, which is used to char- 
acterize the superfluidity of the system, for U — — 4.0i, 
Cd = O.li and filling factor p = 1.2 . The SS melts into 
a normal phase via two successive transitions. First, it 
melts into a phase with zero paring order parameter Aj 
but finite pq at the temperature T ~ 0.2i/fcs, where ks 
is the Boltzmann constant. We will denote this phase 
CDW in the following discussion. One interesting behav- 
ior to be noted in this step of the melting process is that 
as T increases the CDW order parameter keeps grow- 
ing to a maximum until the paring order A decreases to 
zero, which is quite different from the melt ing processes 
of bosonic supersolids investigated in Ref . [17], lis! . The 



physical reason for this behavior is the competition be- 
tween pairing and CDW order in the system. Similar 
phenomena have already been observed in previous stud- 
ies of condensed matter systems, see e.g. Ref. [57|. To 




Figure 3: (Color online). Melting of the supersolid with 
increasing temperature (U = — 4.0t, c<j = Q.lt, p — 1.2). 
Squares and circles correspond to the temperature depen- 
dence of the CDW order pq and the average value of the 
pairing order parameters A, which is used to characterize the 
superfluidity of the system, respectively. 



clarify this point, in the vicinity of the phase transition 
where both order parameters are small, we calculate the 
Ginzburg-Landau (GL) free energy of the system up to 
the 4th order: 



F = A A \A\ 2 + A P \ PQ \ 
+C Ap \A\ 2 \p Q \ 2 . 



B A |A| 4 + i?> Q | 



(6) 



Here {A A , B A , A p , B p , C Ap ] are the GL coefficients 
where a positive C Ap indicates competition between the 
two types of order. A detailed calculation (see appendix 
A) shows that C Ap is indeed positive, thus justifying the 
physical picture given above. Upon further increase of 
the temperature T, the CDW order parameter pq de- 
creases to zero at T ~ Q.29t/kB- 

Moreover, we map out a phase diagram of the system 
with respect to temperature T and the dipolar interaction 
strength Cd for onsite attraction U = — 4M at fixed filling 
p = 1.2 , see Fig. |U When we lower the system tempera- 
ture, for weak dipolar interaction (cd < 0.05i) we observe 
a transition from the normal to the HSF phase. While 
in a very narrow region 0.05t < Cd < 0.07t, we observe 
first a transition between a normal phase (white region) 
and a CDW phase (yellow region), then a transition into 
the SS phase (green region) , and finally a transition into 
the HSF (blue region) which is due to the competition 
between CDW order and pairing as discussed above. At 
larger dipolar interaction strength (cd > 0.07i), we ob- 
serve first a transition between the normal phase and the 
CDW phase, then a transition into the supersolid. As 
can be seen from Fig. 21 in order to observe the super- 
solid for these parameters, Cd should be larger than 0.05<, 
with a maximum critical temperature of about 0.22t/ks 
which is approximately O.ITf . In comparison to the low- 
est temperatures reached in current experiments with po- 
lar molecules in a harmonic trap, which are of order Tp 
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Figure 4: (Color online). Finite temperature phase diagram 
for U — — 4.0t and filling p = 1.2 obtained from calculations 
on a 12 x 12 square lattice. See the text for details. 



[44l | , this means that the temperature still has to be low- 
ered by one order of magnitude in order to observe the 
supersolid. But given major theoretical and experimen- 
tal efforts in lowering the temperature of cold atoms and 
molecules in optical lattices [58l462| . we expect this su- 
persolid phase will be accessible experimentally in the 
near future. 



V. CONCLUSION 

In conclusion, we have shown that when fermionic po- 
lar molecules are loaded into a 2D square lattice, this sys- 
tem will exhibit a supersolid phase which can be observed 
in an experimental setup based on 40 K 87 Rb molecules 
provided the lowest temperature of current experiments 
can be lowered by one order of magnitude. 
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Appendix A: Calculation of GL coefficients 



the Bloch representation, i.e. 
H = 5>(k) - M )cLcka 



kti 



^ V ct c t c ,c. . I 1 V V(q) 

q 



where 



ki k2q 



(Al) 



PqP-q • 



_ xi c k+ q ., 



Ck,, 



k<7 



V(q) = E 



e N 



4 tv q 



|R,f 



(A2) 
(A3) 



and e(k) is Fourier transformation of the hopping matrix. 
The partition function of the system can be written in the 
path- integral form of Z = J T>(cka, Cka)e~ s ^ Ck '" Ck ''} where 
the action of the system is 

S[cka,Ckcr} 

= J~] y^(-io; n + g(k) - /i)cfc a c feo . 



k <y 
U x ^ 

T-pjy 2-^ C ki+q1- C k 2 -ql£k 2 lC kl t 
k±k 2 q 

1 1 

f 2~j3N 



(A4) 



X]^( q ) X C k+q,crCk,aCk' 



-q.p 



Ck' 



kk' aa' 



where w„ = (2n + l)n/(3 is fermionic Matsubara fre- 
quency with f3 — l/(/cgT) being the inverse tempera- 
ture and q, k are three-momentum comprising a vecto- 
rial momentum in 2-dimensional space and a fermionic 
Matsubara frequency. Now by performing two Hubbard- 
Stratonovich transformations 



fci k 2 q 

/^X,x)exp|-E^| X9 | 2 

^X c H c ?^tj + (x 



(A5) 



exp 



Cq-k^Cki Xq 



/ j 2 OjY c k+q,aCk,aCk' -q.cr'Ck' ,a< 

q kk'aa' 



J ^exph^^^y-^q)^ 



(A6) 



In this appendix, we give the detailed calculation of the 
GL coefficients from the microscopic Hamiltonian ([T]) in 



we introduce two auxiliary fields \q an d 4>q which sepa- 
rately relate to the singlet pairing order and charge den- 



6 



sity wave order in the following way: 



(Xq) = -7^(52 c HC q -k t ), 



2i^(q) 

/3N 



E 



(A7) 
(A8) 



After integrating out the fermionic degrees of free- 
dom Cfe CT and Cfc<j, the partition function becomes Z = 
fT>{x,X,<f>)e~ S[x ' x '' l ' ] where 



-trlnCT 1 . 

The matrix Q~ l have the following structure, 



A 



-C7 P 



A 

p-i 



(A9) 



(A10) 



where matrices Gq ~ , A, A. and $ are given by 

^o _1 ) fefe , = (iu n - (e(k) -fJ,))6 k k', Afefe/ = x<fc+fc', 

Afefe' = Xfc+fc', ^/cA:' = -i(0-/c+fc')i and *fefc' = i(0fc-fc')- 

Before presenting a further analysis of the effective ac- 
tion ^[XjX;^] m terms of the specific order parameter 
modes which we are interested in, one thing to be noted 
is that the zero mode of 4> qi i.e. 4>(ofi)i corresponding 
to the particle density of the system, always has a non- 
vanishing contribution to the effective action S[x,Xi 4>] 
independently of the system's parameters (interaction 
strength, temperature, etc). But on the other hand, we 
also note that the effect of </>(o,o) is J us t to renormalize the 
chemical potential /z. Thus in the following analysis we 
simply neglect </>(o,o) m S[x, Xi 4>] by replacing the bare 
chemical potential fj, with a renormalized one /2 which can 
be determined from the particle density of the system. 

To simplify the notation, we further define 



g p-i 








$ 
$ 



,A 



A 
A 



(All) 



(A12) 



We note that 

trlnCT 1 =trln0 o " 1 +t rl n[ 1 + ^oA + ^o*], (A13) 



where the first term just contributes a trivial constant. In 
the vicinity of the phase transition when both the order 
parameter fields p q and A q are small, we can expand the 
second term in (|A13|) in terms of QqA and Go&, i.e., 



trln[l + g A + g $] = J2 



(-1) 



n+l 



-tr(£oA + ao*) 



( Al4 > 

Odd order terms of either A or (j) in the above expan- 
sion vanish since the effective action S'[x,Xi0] preserves 
the symmetry of the unordered phase. Since we are con- 
cerned with finite temperature phase transitions, we ne- 
glect the quantum fluctuations by treating Xq an d 4> q to 
be independent of Matsubara frequency, i.e., we focus on 
the zero Matsubara frequency component of \q an d 4> q - 
Furthermore we assume that X(o,o)i ^(Q.o) an d their con- 
jugate dominate the effective interaction S[x, x, 4>]i where 
Q = (it, it). From the above expansion formula (|A14|1 . 
the Ginzburg-Landau coefficients {Aa, -Ba, A p , B pi Ca p } 
can be calculated explicitly. 

In the following, we calculate the GL coefficient Cap 
which we are most interested in. Expanding S[x, Xi 4 1 ] to 



the 4th order in X(o,o) an d 
X(o,o) a nd 0(q, O ) is given by 



0) 



the cross term between 



2tr 



Go^Go® 



4tr 



GoA 



E 



(£ 2 (k) + uj*)(ioj n - £(k))(iw„ - £(k + Q)) 

|X(O,O)| 2 |0(Q,O)| 2 

(A15) 



(^(k)+ w 2)(^(k+Q)+ W 2) 

C X 4> IX(o,o)| 2 |0(q,o)| 2 



where £(k) = e(k) — /2. It accounts for the compe- 
tition between singlet pairing and CDW order. From 
(|A7|) and (|A8[) we easily obtain (x(o,o)) = W\A an d 
(cj) q ) = 2iV(q)p Q which gives C Ap = C xl pW 2 (Q)U 2 
after identifying the effective action S[x, Xj</>] with the 
Ginzburg-Landau free energy. 

Now we calculate C x< p. The summation over Matsub- 
ara frequency can be performed explicitly: 
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E 



(e 2 (k)+^ )( ^„_ e(k))(lWn _ e(k + Q)) - (£2 (k)+w 2 )(e 2 (k+Q)+w 2 ) 

4 2 



y 2ni'T eP* + l 



1 



(z-£(k)) 2 (z + £(k))(z-£(k + Q)) (z 2 -e(k))(z 2 -e(k + Q)). 



+ 



f (k) + £(k + Q) V (1 + e«( k +Q))(£(k + Q) - £(k)) 2 (1 + e' 3 «( k ))^ 2 (k) 
g(k + Q) - 3g(k) + e^( k )(g(k + Q) + C(k)(-3 + 2/3(g(k + Q) - g(k)))) 
(1 + e /3«( k )) 2 ^(k)(^(k + Q) - £(k)) 2 
1 AanhiOs+Qi tanh«M^ 



C 2 (k)-e 2 (k + Q) ^ £(k + Q) 



e(k) 



(A16) 



200 
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Figure 5: Temperature dependence of C x< f, at filling p — 1.2 . 



The remaining summation over momentum k can be eas- 
ily replaced by an integral over the density of states 
Pdos(^) = Sk^( e — e (k)) an d evaluated numerically. 
As we can see from Fig. O which shows the tempera- 
ture dependence of C x ^ when 0.1 < ksT/t < 0.3 at fill- 
ing p = 1.2, C xt /) is positive in this temperature region. 
Therefore, since the GL coefficient Ca p has the same sign 
as C x< p , we now clearly see that singlet pairing and CDW 
order indeed compete with each other in the parameter 
region under investigation. 
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